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Abstract 

Discontinuous dynamical systems with grazing solutions are discussed. The group property, continu¬ 
ation of solutions, continuity and smoothness of motions are thoroughly analyzed. A variational system 
around a grazing solution which depends on near solutions is constructed. Orbital stability of grazing 
cycles is examined by linearization. Small parameter method is extended for analysis of neighborhoods 
of grazing orbits, and grazing bifurcation of cycles is observed in an example. Linearization around an 
equilibrium grazing point is discussed. The mathematical background of the study relies on the theory 
of discontinuous dynamical systems [T]. Our approach is analogous to that one of the continuous dy¬ 
namics analysis and results can be extended on functional differential, partial differential equations and 
others. Appropriate illustrations with grazing limit cycles and bifurcations are depicted to support the 
theoretical results. 

Keywords: Discontinuous dynamical systems; Grazing points and orbits; Axial and non-axial grazing; 
Variational system; Orbital stability; Small parameter; Bifurcation of cycles; Impact mechanisms 

1 Introduction 

Vibro-impacting systems provide examples of non-linear dynamical systems, exhibiting new levels of 
complicated dynamics due to their non-smoothness. Grazing phenomenon is one of the attractive features 
for these dynamics Ca-ES], M-M- There are two approaches for the definition of grazing in literature. 
One is presented in the studies of Bernardo, Budd and Ghampneys EIH], Bernardo and Hogan [a, and 
Luo [lS]-[in|- these studies, it is asserted that grazing occurs when a trajectory hits the surface of 
discontinuity tangentially. In [SSJ-ESli Nordmark defines grazing as the approach of the velocity to zero 
in the neighborhood of the surface of discontinuity which is the case of the studies conducted in m-m, 
|26|-|29j. Our comprehension of grazing in this paper is close to that one in [H],ia,12S]. In the paper m, 
grazing is considered as a bounding case which separates regions of quite different dynamic behaviors. It 
is understood that the system trajectory makes tangential contact with an event triggering hypersurface. 
The shooting, continuation and optimization methods are developed and illustrated for both transient 
and grazing phenomena. It is exemplified by utilizing power electronics and robotics. In m, the 
grazing periodic orbit and its linearization are obtained by means of a numerical continuation method 
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for hybrid systems. Applying this, the normal-form coefficients are evaluated, which in this case imply 
the occurrences a jump to chaos and period-adding cascade. The necessary and sufficient conditions 
of the general discontinuous boundary are expounded in |28| . In mi: by means of non-stick mapping, 
the necessary and sufficient conditions for the grazing in periodically forced linear oscillator with dry 
friction are obtained. By constructing special maps such as zero time discontinuity mapping m-m 
and Nordmark map |36|-|39|. the existence of periodic solution and their stability were investigated in 
mechanical systems. 

In this paper, we model the dynamics with grazing impacts by utilizing differential equations with 
impulses at variable moments and applying the methods of m-m- As a consequence of such methods, 
the role of the mappings |36|-|39j is diminished. One can observe that a trajectory at a grazing point 
may have tangency to the surface of discontinuity, which is parallel to one or several coordinate axises. 
Particularly it means the velocity approaches to zero [3S]-[3ni- Then, we will say about the axial grazing. 
Otherwise, grazing is non-axial. This research contains the analysis of both axial and non-axial grazing. 

In [221, it is observed through simulations and experiments that the coefficient of restitution depends 
on the impact velocity of the particle by considering both the viscoelastic and the plastic deformations 
of particles occurring at low and high velocities, respectively. It has been proposed in m that at low 
impact velocities and for most materials with linear elastic range, the coefficient of restitution is of the 
form R{v) = I — av, where v is the velocity before collision and a is a constant. Also, for low impact 
velocities the restitution law can be considered quadratic m- This is why, we will use non-constant 
restitution coefficients in models with impacts in this study. 

For investigation of autonomous differential equations, it is convenient to utilize properties of dynam¬ 
ical systems. They are the group property, continuation of solutions in both time directions, continuity 
and differentiability in parameters. The studies of discontinuous dynamical systems with transversal 
intersections of orbits and surfaces, B— smooth discontinuous flows, can be found in mm- In this 
research, the dynamics is approved for systems with grazing orbits. Moreover, the definitions of orbital 
stability and asymptotic phase are adapted to grazing cycles. The orbital stability theorem is proved, 
which can not be underestimated for theory of impact mechanisms. 

The remaining part of the paper is organized as follows. In the next section, we will introduce 
necessary notations, definitions and theorems to specify discontinuous dynamical systems. In Section 
3, it is shown how the dynamics can be linearized around grazing orbits. In Section 4, the theorem of 
orbital stability is adapted for grazing cycles. In Section 5, the small parameter analysis is applied near 
grazing orbits and bifurcation of cycles is observed. In last three sections, examples are presented to 
actualize the theoretical results numerically and analytically. Finally, Conclusion covers a summary of 
our study. 
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2 Discontinuous dynamical systems 


Let K., N and 1 be the sets of all real numbers, natural numbers and integers, respectively. Consider the 
set D G M" such that D = UDi, where Di, i = 1,2,... ,k, components of D, are disjoint open connected 
subsets of R". To describe the surface of discontinuity, we present a two times continuously differentiable 
function $ : —>■ R". The set can be defined as T = $“^(0) and is a closed subset of D, where D is the 

closure of D. Denote 9r as the boundary of F. One can easily see that T = where F^ are parts 

of the surface of discontinuity in the components of D. Denote F = J(T), $(x) = $(J“^(a::)). Denote an 
r— neighborhood of D in R" for a fixed r > 0 as D^. Let F’’ be the r— neighborhood of F in R”, for a 
fixed r > 0 and define functions J : F’' —?► O’" and J : f , such that, J(F), J(f) C D. Assume that 
a function f{x) : —>■ R" is continuously differentiable in . Set the gradient vector of $ as 

The following definitions will be utilized in the remaining part of the paper. Let x(t—) be the left 
limit position of the trajectory and x{t+) be the right limit of the position of the trajectory at the 
moment t. Define Ax{t) := x{t+) — x(t—) as the jump operator for a function x{t) such that x{t) G F 
and t is a moment of discontinuity (discontinuity moment). In other words, the discontinuity moment t 
is the moment when the trajectory meets the surface of discontinuity F. The function I{x) will be used 
in the following part of the paper which is defined as I(x) := J{x) — x, for a: € F. 

The following assumptions are needed throughout this paper. 


(Cl) 

(C2) 

(C3) 

(C4) 

(C5) 

(C6) 

(C7) 


V$(x) ^ 0 for all a: £ F, 

J £ C'i(F'') and det ^ 0, for all a: £ F’' \ ^F, 


Ffif c aFnar, 

(V$(x), /(a:)) ^ 0 if X £ F \ aF, 

(vl>(x), /(x)) ^ 0 if X £ f \ af, 

J(x) = X for all X £ 3F, 

J(x) = X for all X £ 9f. 


One can verify that f = {x £ D|cE>(x) = 0} and J(x) x on f since of {C2). Condition (Cl) implies 
that for every xq £ F, there exist a number j and a function (f)xQ (xi,..., Xj_i, Xj+i,..., x„) such that F 
is the graph of the function Xj = ■. ■ ,Xj_i,Xj+i,... ,x„) in a neighborhood of xq. Same is true 

for every xq £ F. Moreover, V$(x) ^ 0, for all x £ F, can be verified by using the condition (C2). The 
conditions (C2), (C6), (C7), imply that the equality J{x) = x, is true for all x £ clF. 

Let .e/ be an interval in Z. We say that the strictly ordered set 9 = {9i},i G £/, is a B—sequence 
[T] if one of the following alternatives holds: ii) 0 = 0, {ii) 0 is a nonempty and finite set, (Hi) 9 is an 
infinite set such that \9i\ -G 00 as f > 00 . In what follows, 0 is assumed to be a i3—sequence . 
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The main object of our discussion is the following system, 


x' = f{x), 
Aa;|a;Gr = d(a;). 


( 2 . 1 ) 


In order to define a solution of (EB, we need the following function and spaces. 

A function : R —>■ R", n G N, d is a il—sequence, is from the set PC'(R, 0) if it : (i) is left 

continuous, (ii) is continuous, except, possibly, points of 6, where it has discontinuities of the first kind. 

A function is from the set PC^(R^0) if (j)(t), (j)'(t) G PC(JS.,9), where the derivative at points of 
0 is assumed to be the left derivative. If (j){t) is a solution of (12.11) . then it is required that it belongs to 
PC\R,0) p. 

We say that x{t) : —>■ R", C R, is a solution of (12.11) on if there exists an extension x{t) of 
the function on R such that x(t) G ^(^^(R,0), the equality x'(t) = f{x(t)), t G is true if x{t) ^ T, 
x{0i+) = J{x{0i)) for x{0i) G r and x{0i+) G T, G J^. If 0i is a discontinuity moment of x{t), then 

x{0i) G r, for 0i > 0 and x{0i) G f, for 0i < 0. If x{0i) G ST or x{0i) G (?f , then x{0i) is a point of 

discontinuity with zero jump. 

Definition 2.1 A point x* from idT or idT is a grazing point of system (12.11) if (V<i)(a;*), f{x*)) = 0 or 
(V<P(x*), f(x*)} = 0, respectively. If at least one of coordinates of V^ix*) is zero then the grazing is 
axial, otherwise it is non-axial. 

Definition 2.2 A n orbit j(x*) = {a;(t, 0, a;*)|x* G D, t G R} of \2.1\) is grazing if there exists at least 
one grazing point on the orbit. 

Consider a solution x{t) : R —>■ R" and {0i} be the moments of the discontinuity, they are the moments 
where solution x{t) intersects F as time increases and the moments when the solution it intersects F as 
time decreases. 

A solution x{t) = x{t, 0, xq), xq G D of (12.11) locally exists and is unique if the conditions (Cl) — (C3) 
are valid [T]. 

In what follows, let || • || be the Euclidean norm, that is for a vector x = (xi,X 2 ,... ,x„) in R", the 
norm is equal to y^x^ + x| + ... + x^. 

The following condition for m guarantees that any set of discontinuity moments of the system 
constitutes a B— sequence and we call the condition B— seguence condition. 

(C8) sup£, ||/(x)|| < + 00 , and inf,^^gp(xo, y(C, 0, Xq)) > 0. 

In [Tj, some other B— sequence conditions are provided. 

We will request for discontinuous dynamical systems that any sequence of discontinuity moments to 
be a P— sequence. 
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Let us set the system 


y' = fiy) 


( 2 . 2 ) 


for the possible usage in the remaining part of the paper. 

Consider a solution y{t,0,xo), xq € L, of (12.211 . Denote the first meeting point of the solution with 
the surface L, provided the point exists, by ?/((C, 0,xo). The following conditions are sufficient for the 
continuation property. 

(C9) (a) Every solution y{t,0,xo), xq € D, of (12.21) is continuable to either oo or T as time increases, 
(b) Every solution y{t, 0,xq), Xq € D, of (j2.2p is continuable to either — oo or E as time decreases. 

To verify the continuation of the solutions of (EU, the following theorems can be applied. 

Theorem 2.1 m Assume that conditions (CS) and ((79) are valid. Then, every solution x{t) = 
x{t,0,xo), xo € D of \2.1\) is continuable on R. 

Now, we will present a condition which is sufficient for the group property. 

(CIO) For all xq G D, the solution y{t, 0, xq) of (12.21) does not intersect F before it meets the surface T as 
time increases. 

In other words, for each Xq G D and a positive number s such that y(s, 0, xq) G T, there exists a number 
r, 0 < r < s, such that y{r, 0, Xq) G F. 

It is easy to verify that the condition ((710) is equivalent to the assertion that for all xq G D, the 
solution y{t,0,xo) of (12.21) does not intersect F before it meets the surface F as time decreases. In other 
words, for each xq G D and a negative number s such that y(s, 0,xo) G F, there exists a number r, 
s < r < 0, such that y{r, 0, xq) G F. 

Theorem 2.2 fThe group property) Assume that conditions (Cl)-(CIO) are valid. Then, 
x(t 2 ,0, x(ti, 0, xo)) = x{t 2 +ti,0, xo), for all ti,t 2 G K. 

Proof. Denote by ^(t) = x{t + f), for a fixed t G R. It can be verified that the sequence {6i — t} is a set 
of discontinuity moments of f{t) and the function is a solution of (12.11) [T]. The next step is to show that 
the following equality x(—t, 0, x(t, 0, Xq)) = xq, holds for all xq G D and t G R. Consider the case t > 0. 
If the set of discontinuity moments {0i} is empty, the proof is same with that for continuous dynamical 
systems |23| . Because of the condition ((72), which corresponds to invertibility of the jump function J, 
the equality x{6i,0,x{9i+)) = x{6i), holds for all i G £/. Assuming that 6(_i < 0 < 0i, we should verify 
x(—01,0, x(6(i, 0, Xo)) = Xq. Denote by x(t) = x(t, 0, x(0i)). The point x(0i) lies on the discontinuity 
surface F. By condition ((73) the solution x(t) is a trajectory of y' = /(y) for decreasing t. Condition 
((710), part (a), implies that the trajectory x{t) cannot meet with F if t > —9i as time decreases. That 
is, x(—0i) = Xo as the dynamics is continuous. The proof for t < 0 can be done in a similar way. □ 
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Remark 2.1 For the application of the results, it is possible to take the initial moment as to = 0, without 
being the discontinuity moment since of the group property. Then Xq ^ F U F. 

Denote by [a,b], a,b £ M, the interval [a,b], whenever a < b and [b,a\, otherwise. Let a;i(t) £ 
PC{M.+ ,6^), 6^ = {0}}, and X 2 {t) £ PC{M.+ ,6‘^), 0^ = {Of}, be two different solutions of (12.11) . 

Definition 2.3 The solution X 2 {t) is in the e—neighborhood of xi{t) on the interval if 

• the sets 6^ and 6^ have same number of elements in 

• ~ (^i\ < e for all 9] £ 

• the inequality ||xi(t) — X 2 (t)\ \ < e is valid for all t, which satisfy t ^ \ Ugi 

The topology defined with the help of e— neighborhoods is called the B-topology. It can be apparently 
seen that it is Hausdorff and it can be considered also if two solutions Xi{t) and X 2 (t) are defined on a 
semi-axis or on the entire real axis. 

Definition 2.4 The solution XQ{t) = x{t,0,xo), t € R, xq £ D, of h2.1\) B-continuously depends on xq 
for increasing t if there corresponds a positive number 5 to any positive e and a finite interval [0, b] , b>0 
such that any other solution x(t) = x{t,0,x) of h2.1\) lies in e—neighborhood of xo(t) on [ 0 , 6 ] whenever 
x £ B{xo,S). Similarly, the solution Xo(t) of i2.1\) B-continuously depends on Xq for decreasing t if there 
corresponds a positive number 6 to any positive e and a finite interval [a, 0], a < 0 such that any other 
solution x(t) = x(t,0,x) of i2.1\) lies in e—neighborhood of x^it) on [a,0] whenever x £ B(xo,S). The 
solution Xo{t) of i2.1\) B-continuously depends on xq if it continuously depends on the initial value, Xq, 
for both increasing and decreasing t. 

If conditions (CI)-(C7) hold, then each solution xo{t) : K. —>• M", xo{t) = x{t,0,xo), of ( 12 . 11 ) continu¬ 
ously depends on xq [T]. 

2.1 B-equivalence to a system with fixed moments of impulses 

In order to facilitate the analysis of the system with variable moments of impulses m, a B-equivalent 
system [1] to the system with variable moments of impulses will be utilized in our study. Below, we will 
construct the B-equivalent system. 

Let x{t) = x{t,0,xo + Ax) be a solution of system ( 12 . 11 ) neighbor to xo{t) with small ||Aa:||. If the 
point xo(Oi) is a (/?)— or ( 7 )— type point, then it is a boundary point. For this reason, there exist two 
different possibilities for the near solution x(t) with respect to the surface of discontinuity. They are: 

(iVl) The solution x(t) intersects the surface of discontinuity, F, at a moment near to 0i, 

(N2) The solution x(t) does not intersect F, in a small time interval centered at 0^. 
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Consider a solution xo(t) : J" R", C M., oi (12.11) . Assume that all discontinuity points 9i, 
i £ £/ are interior points of . There exists a positive number r, such that r-neighborhoods of Di[r) of 
{Oi,xo{6i)) do not intersect each other. Consider r is sufficiently small and so that every solution of (12.21) 
which satisfies condition (A^l) and starts in Di{r) intersects F in G'i(r) as t increases or decreases. Fix 
i £ s/ and let ^{t) = x{t,9i,x), {9i,x) £ Di{r), be a solution of (12.21) . Tj = Ti{x) the meeting time of ^{t) 
with r and = x{t,Ti,^{Ti) + Ji^in))) another solution of (12.21) . Denoting by Wi{x) = '4’{9i) — x, 
one can find that it is equal to 

W^{x)=[ f{^is))ds + J{x+ [ f{^{s))ds)+[ /(V'(s))ds (2.3) 

J0i JOi Jn 

and maps an intersection of the plane t = 9i with Di{r) into the plane t = 9i. 

Let us present the following system of differential equations with impulses at fixed moments, whose 
impulse moments, {9i}, i £ si^ are the moments of discontinuity of Xo{t)^ 


y' = fiy). 


Ay\t^g. = Wi{y{9,)). 


(2.4) 


The function / is the same as the function in system and the maps Wi, i £ si, are defined by 
equation (12.31) . if ^(t) = x{t,9i,x) does not intersect T near 9i then we take Wi{x) = 0. 

Let us introduce the sets = {{t,x)\t G I, ||x — xo(t)|l < r}, and Di{r), i £ si, closure of an 
r— neighborhood of the point {9i,xo{9i+)). Write D’’ = F). U (Uigj^^Di(r)) U (Uigj^^Di(r)). Take r > 0 
sufficiently small so that D’’ C R x D. Denote by D{h) an h-neighborhood of a;o(0). Assume that 
conditions (Cl) — (CIO) hold. Then systems (12.11) and (12.41) are B-equivalent in for a sufficiently 
small r [T]. That is, if there exists h > 0, such that: 


1. for every solution y{t) of (12.41) such that y{0) £ D{h), the integral curve of y{t) belongs to and 
there exists a solution x{t) = x{t,0,y{0)) of ( 12 . 11 ) which satisfies 


x{t) = yit), t £ [a,&]\U™_fc {Ti,9,], 


(2.5) 


where are moments of discontinuity of x(t). One should precise that we assume 
satisfies {N2). Particularly, 


x{9i) 


y{9i), if 9i < Ti, 

< 

2 /( 0 )*'), otherwise, 

► 


V{Ti) 


x{Ti), if 9i > Ti, 
x{t^), otherwise. 


9i, if x{t) 


( 2 . 6 ) 


2. Conversely, if (|2.4D has a solution y{t) = y(t, 0, t/(0)), 2 /( 0 ) £ D(h), then there exists a solution 
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x(t) = x{t,0,y{0)) of (12.11) which has an integral curve in , and (12.51) holds. 


A solution xo{t) satisfies (12.11) and (12.41) simultaneously. 

Consider a solution xo{t) : R —R”, xo{t) = x{t,0,xo), xq G D with discontinuity moments {Oi}. 
Fix a discontinuity moment 9i. At this discontinuity moment, the trajectory may be on F and F. All 
possibilities of discontinuity moment should be analyzed. For this reason, we should investigate the 
following six cases: 

(a) XQ{9i) e F \ 9F, (o') xo{9i) G f \ 9f, 

(^) xo{9,) G aF & (V$(xo(0.)),/(^o(0.))) ^ 0, (/?') x^{9,) G af & (V$(a:o(0z)),/(xo(0*))) 0, 

( 7 ) xo(0O G aF & (V$(xo(dO),/(^o(0*))) =0, (Y) xo(Y) e af & (V$(xo(Y)),/Yo(Y))) =0. 

If a discontinuity point xo{9i) satisfy the case (a), {{a')) the case (/3), ((/?')) and the case ( 7 ), (( 7 ')) 
we will call it an (a)— type point, a (/3)— type point and a ( 7 )— type point, respectively. 

Besides, we present the following definition which is compliant with Definition 12.21 

Definition 2.5 If there exists a discontinuity moment, 9i, i G for which one of the cases ( 7 ) or ( 7 ') 
is valid, then the solution xo(t) = x(t, 0,a:o), xq G R" of i2.1\} is called a grazing solution and t = 9i is 
called a grazing moment. 

Next, we consider the differentiability properties of grazing solutions. The theory for the smoothness 
of discontinuous dynamical systems’ solutions without grazing phenomenon is provided in [1]. 

Denote hy x{t),j = 1, 2,..., n, a solution of (12.41) such that a;(0) = xq + Ax, Ax = (fy, fy, ■ ■., fy), 
and let pi be the moments of discontinuity of x(t). 

The following conditions are required in what follows. 

(A) For all t G [0, b]\ {rji, Y], the following equality is satisfied 

n 

x{t) - xo{t) = ^ Ui{t)f,^ + 0(11 Ax||), (2.7) 

where Ui(t) G PC{[0,b],9). 

(B) There exist constants I'ijjj G £/, such that 

n 

-^3 = + 0(11 Ax||); (2.8) 

i=l 

(O) The discontinuity moment pj of the near solution approaches to the discontinuity moment 9j,j G 
of grazing one as tends to zero. 
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The solution x{t) has a linerization with respect to solution xo{t) if the condition (A) is valid and, 
moreover, if the point Xo(ffi) is of (a)— or (/3)— type, then the condition (B) is fulfilled. For the case 
xo(0i) is of ( 7 )— type the condition (C) is true. 

The solution xo(t) is iF—differentiable with respect to the initial value xq on [0, h] if for each solution 
x(t) with sufficiently small Aa; the linearization exists. The functions Ui{t) and Vij depend on Ax and 
uniformly bounded on a neighborhood of xq. 

It is easy to see that the differentiability implies S—continuous dependence on solutions to initial 
data. 

Define the map C{t,x) as C{t,x) = x{t,0,x), for x G D. 

A AT-smooth discontinuous flow is a map ({t, x) : R x D —>■ D, which satisfies the following properties: 

(I) The group property: 

(i) C(Oj^) ■ D ^ D is the identity; 

(ii) ({t, ({s, x)) = Q{t + s, x) is valid for all t, s S R and x G D. 

(II) C(^)^) G PC'^(R) for each fixed x G D. 

(Ill) ({t, x) is AT-differentiable in x G D on [a, 6] C R for each a, b such that the discontinuity points of 
C(t,x) are interior points of [a, b]. 

In [1], it was proved that if the conditions of Theorem 12.11 and (CI)-(CIO) are fulfilled, then system 
(EH) defines a D-smooth discontinuous flow [1] if there is no grazing points for the dynamics. It is easy 
to observe that the i?-smooth discontinuous flow is a subcase of the AT-smooth discontinuous flow. In 
the next section, we will construct a variational system for EH in the neighborhood of grazing orbits. 
That is, we will assume that some of the discontinuity points are ( 7 )— type points. Linearization around 
a solution and its stability will be taken into account. Thus, analysis of the discontinuous dynamical 
systems with grazing points will be completed. 

3 Linearization around grazing orbits and discontinuous dynam¬ 
ics 

The object of this section is to verify K— differentiability of the grazing solution. Consider a grazing 
solution XQ{t) = x{t,0,xo), Xq G D, of (12.11) . We will demonstrate that one can write the variational 
system for the solution xo{t) as follows: 


u' = A{t)u^ 
Au\t=g, = 


(3.9) 
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where the matrix A{t) € R”^" of the form A(t) = The matrices i = 1,... ,n, will be 

defined in the remaining part of the paper. The matrix Bi is bivalued if 0i is a grazing moment or of 
(/3)-type. 

The right hand side of the second equation in (13.911 will be described in the remaining part of the 
paper for each type of the points. As the linearization at a point of discontinuity, we comprehend the 
second eqnation in (13.911 . 


3.1 Linearization at (a)— type points 

Discontinuity points of (a) and {a') types are discussed in [I]. In this subsection, we will outline the 
results of the book. 

Assume that x{9i) is an (a)—type point. It is clear that the B— equivalent system (12.411 can be 
applied in the analysis. The functions Ti(x) and Wi{x), are described in Subsection 12.II Differentiating 
<I>(a;(Ti(a;))) = 0 , we have 

dn{xo{ei)) 

dxj ^x{xo{9i))f{xQ{9i))' 

Then, considering (12.3L we get the following equation, 


dW,{x^{9i)) 


dxi 


oj 


= {f{xo{9i)) - /(xo( 6 '*) + J{xo{9,)))) 


dxi 


dxi 


Oj 


dx 
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(3.11) 


where tj = (0 ,..., 1 ,..., 0). 

3 _ 

The matrix Bi G R”^” in equation (13.911 is defined as Bi = Wix, where Wix is the nxn matrix of the 
form Wix = dWi{xo{ei)) ^ ^ vector-components , j = 1 ,..., n, 

evaluated by (13.1111 . Moreover, the components of the gradient Vr^ have to be evaluated by formula 

dnni). 


3.2 Linearization at (/?)—type points 


In what follows, denote nxn zero matrix by On- In the light of the possibilities (A^l) and {N2), the 
matrix Bi in ( 12.111 can be expressed as follows: 


B^= < 


On, if 


W,x, if 


(7V1) is valid, 
(7V2) is valid. 


(3.12) 


where Wix is evaluated by formula (13.1111 and Vt(x) evaluated by formula (13.1011 . 

The differentiability properties for the cases (a') and (/3') can be investigated similarly. 
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3.3 Linearization at a grazing point 

Fix a discontinuity moment 6 i and assume that one of the cases ( 7 ) or ( 7 ') is satisfied. We will investigate 
the case ( 7 ). The case ( 7 ') can be considered in a similar way. 

Considering condition (Cl) with the formula (13.101) . it is easy to see that one coordinate of it is 
infinity at a grazing point. This gives arise singularity in the system, which makes the analysis harder 
and the dynamics complex. Through the formula (13.101) . one can see that the singularity is just caused 
by the position of the vector field with respect to the surface of discontinuity and the impact component 
of the dynamical system does not participate in the appearance of the singularity. To handle with the 
singularity, we will rely on the following conditions. 

(Al) A grazing point is isolated. That is, there is a neighborhood of the point with no other grazing 
points. 

{A2) The map Wi{x) in (12.31) is differentiable at the grazing point x = Xo{9i). 

(A3) The function Ti{x) does not exceed a positive number less than 0^+1 — 9i near a grazing point, 
XQ{9i), on a set of points which satisfy condition (iVl). 

In the present paper, we analyze the case, when the impact functions neutralize the singularity 
caused by transversality. That is, the triad: impact law, the surface of discontinuity and the vector 
field is specially chosen, such that condition (A2) is valid. Presumably, if there is no of this type of 
suppressing, complex dynamics near the grazing motions may appear [3 [Ml ISa [37]. In the examples 
stated in the remaining part of the paper, one can see the verification of (A2), in details. 

Let us prove the following assertion. 

Lemma 3.1 If conditions (Cl), (C4), (C 6 ), (C 8 ) and (A3) hold. Then, Ti{x) is continuous near a 
grazing point xo{9i), on a set of points, which satisfy condition (iVl). 

Proof. Let xo(9i) be a grazing point. If x is not a point from the orbit of the grazing solution, 
the continuity of Ti(x) at the point x = x can be proven using similar technique presented in [T]. Now, 
the continuity at xo(9i) is taken into account. On the contrary, assume that Ti(x) is not continuous 
at the point x = xo(9i). Then, there exists a positive number eo and a sequence {xnjnez such that 
Ti{xn) > eo whenever Xn —>• Xci{9i), as n —> 00 . Moreover, from condition (A3), one can assert that there 
exists a subsequence Ti{xnk) which converges to a number eg < tq < 0^+1 — 9i. Without loss of generality, 
assume that the subsequence converges the point where the sequence converges. Since of the 

continuity of solutions in initial value, x{Ti{xn),0,Xn) approaches to x(To,0,x(9i)). But x(Ti(xn),0,Xn) 
is on the surface of discontinuity L, x(ro, 0, a:o(0i)) ^ L. This contradicts with the closeness of the surface 
of discontinuity T. The continuity at other points of the grazing orbit is valid by the group property. □ 
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Since of S—equivalence of systems m and (1^ . we will consider linearization around as 

solution of the system (12.41) . consequently, only formula (12.71) will be needed. Finally, the linearization 
matrix for the grazing point also has to be defined by the formula (|3.12D , where Wix exists by condition 
(A2). 

In what follows, we will consider only grazing motions such that condition {A2) holds. Consequently, 
the continuous dependence on initial data is valid. More precisely, B— continuous dependence on initial 
data is true. Now, if conditions (Cl) — (CIO) and (Al), {A2) are assumed, the system (12.11) defines a 
K— smooth discontinuous flow for dynamics with grazing points. 


3.4 Linearization around a grazing periodic solution 


Let 4('(t) : IR. —>■ U be a periodic solution of (12.11) with period a; > 0 and 6 i, i £ Z, are the points of 
discontinuity which satisfy (uj,p)— property, i.e. + a;,p is a natural number. 

Let us fix a solution x{t) = a:(t, 0, 4'(0) + Aa:) and assume that linearization of 4('(t) with respect to 
x{t) exists and is of the form 


u' = A{t)u^ 

= BiU. 


(3.13) 


The matrix Bi is determined by (13.121) . It is known that A{t + w) = A{t), f £ R. But, the sequence Bi 
may not be periodic in general, since of (13.121) . This makes the analysis of the neighborhood of 4'(t) 
difficult. For this reason, we suggest the following condition. 


(A4) For each sufficiently small Ax £ R", the variational system (13.131) satisfies Bi+p = Bi, i & There 
exist a finite number m < 2}, where I is the number of points of (/?)— or ( 7 )— type in the interval 
[ 0 , 0 ;], of the periodic sequences Bi. 

The assumption (44) is valid for many low dimensional models of mechanics and those which can be 
decomposed into low dimensional subsystems. To distinguish periodic sequences Bi in the assumption 
(44), we will apply the notation Bi = i £ Z and j = 1, 2,..., m. 

If the condition (44) is not fulfilled, then complex dynamics near a periodic motion may appear. 
This case can be investigated either by methods developed through mappings applications [nin] or it 
requests additional development of our present results. 

In the next example, we will demonstrate that the system constitutes K— smooth discontinuous flow 
although it has grazing points in the phase space. 

Example 3.1 (K-smooth discontinuous flow with grazing points). Consider an impact model 


y'l = y2, 


1/2 = - 2/1 + 0 . 0012 / 2 , 


(3.14a) 
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A2/2|yeri = -yi - Rivl 


( 3 . 14 b) 


Ay2|yer2 = -(1 + -^ 2 ) 2 / 2 , 

with the domain D = Ri = exp(—O.OOOS tt) and R 2 = 0.9. In the paper (16]/ . it is stated that the 
coefficient of restitution for low velocity impact still remains as an open problem. In the study m by 
eonsidering Kelvin- Voigt model for the elastic impact, we derived quadratic terms of the velocity in the 
impact law. This arguments make the quadratic term for the impulse equation (I3.14bl) reasonable. 

Let us describe the set of discontinuity curves by T = riUr 2 . The eomponents Fi and r 2 are 
intervals of the vertical lines yi = exp(0.000257r) and yi = 0, respectively and they will be precised next. 
Fix a point P = (0, j/ 2 ) € D, with 1/2 > 1. Let y(t,0,P) be a solution of and it meets with the 

vertical line xi = exp(0.000257r), X 2 > 0 at the point P 2 = (exp(0.000257r), 2 / 2 ( 6 *!, 0, P)), where 9i is the 
meeting moment with the line. Consider the point Q 2 = (exp(0.000257r), —i?i 2 / 2 ( 0 i, 0, P 2 )^) OLnd denote 
Qi = (0, 2 / 2 (^ 2 ,0, ( 32 )), where 02 is the moment of meeting of the solution y(t,Q,Q 2 ) with the vertical line 
Xi =0, X 2 < 0. We shall need also the point Pi = (0, —P 2 y 2 (^ 2 , 0, (32))- Finally, we obtain the region 
G in yellow and blue between the vertical lines and graphs of the solutions in Figure [IJ The region G 
eontains diseontinuous trajectories and outside of this region all trajectories are continuous. Moreover, 
both region G and its complement are invariant. 



Figure 1 : The region G for system ( 13 . 141 ) is depicted in details. The curves of discontinuity T = Ti U r2 
and r = Ti U r2 are drawn as vertical lines in red and green, respectively and the grazing orbit in 
magenta. 


Define Ti = {{yi,y2)\ yi = exp( 0 . 000257 r), 0 < 2/2 < 2/2(6*!, 0 , ( 0 ,2/2))}, and r2 = {(2/!, 2/2)! yi = 

0 , 2/2(6*2, 0 , —Pi2/2((0!, 0 , ( 0 , ^2)))^) < 2/2 < 0 }. The boundary of the curve, T = Ti U r2, has of four 
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points, they are 


dr = {( 0 , 0 ), (exp( 0 . 00025 ^), 0 ), (exp( 0 . 00025 ^), y 2 { 0 i, 0 , ( 0 , ^2))), ( 0 ,2/2(02, 0 , -i?i2/2(0i, 0 , ( 0 , ya)))}- 

In the following part of the example, we will show that two of them, y* = ( 2 /*, 2 / 2 ) = (exp( 0 . 000257 r), 0 ) 
and the origin, ( 0 , 0 ) are grazing points. Moreover, it can be easily validated that other two points are of 
13-type. 

Issuing from system \,‘I.14\ ), the curve of diseontinuity F consists of two components F 1 and F 2 . The 
components are the following sets 

fi = {(2/1,2/2)1 2/1 = exp( 0 . 0002571 ), -i?i2/2(0i, 0 , (0,^2))^ < 2/2 < 0 } 


and 

^2 = {(2/1,2/2)1 2/1 = 0, 0 < --^ 22 / 2 ( 02 , 0 ,( 32 )}- 
One can verify that the function 

{ exp(0.00052) fsin(i),cos(2)V ift£ [ 0 , 7 r), 

^ ^ (3.15) 

(0,1), ift = TT, 

is a discontinuous periodic solution of with period w = tt, whose discontinuity points ( 0 , 1 ) and 

(0, — exp(0.00057r)) belong to F and F, respectively. The expression 

(V$((exp(0.000257r), 0)), /((exp(0.000257r), 0))) = ((1, 0), (0, - exp(-0.000257r))) = 0 


verifies that y* is a ( 7 )— type point, i.e. a grazing point of the solution (2). It is easily seen that the 
grazing is axial. Now, we can assert that the periodic solution SS.15\) is a grazing solution in the sense 
of DeRnition \2.,^ Its simulation is depicted in Figure\^ 



Figure 2: The grazing orbit of system (I3.14l> . 


Since the complement of G is invariant in both directions and consists of continuous trajectories of 
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the linear system I3.14a\ ), one can easily conclude that the complement is a continuous dynamical system 
Thus, to verify the dynamics for the whole system, one need to analyze it in the region G. This set 
is hounded, consequently for solutions in it conditions (C 8 ) and (CQ) are fulfilled and by Theorem \2.1[ 
they admit B—sequences and continuation property. 

Consider a function C( 2 / 2 ) : [ 2 / 2 ( 02 , 0, ( 52 ), 2 / 2 ( 01 ,0, P)] —[ 2 / 2 ( 02 ,0, ( 52 ), 2 / 2 ( 0 i, 0, P)] such that it is 
continuously differentiable, satisfies ^( 2 / 2 ) = ~P 22/2 in a neighborhood 0 / 1/2 = 0 and is the identity at 
the boundary points, i.e. ^( 2 / 2 ( 01 ,0, P)) = 2 / 2 ( 01 ,0,P) and C( 2 / 2 ( 0 i, 0, P)) = 2 / 2 ( 01 , 0 , P). It is easily seen 
that such function exists. On the basis of this discussion, let us introduce the following system, 

y'l = 2/2, 

2/2 = “2/1 + O.OOI2/2, ( 3 . 16 ) 

A 2 / 2 |ygr = C( 2 / 2 ) - 2 / 2 - 

It is apparent that system S3.16\) is equivalent to near the orbit of periodic solution '!'(/). That 

is, they have the same trajectories there. 

Specifying a for iS.16\} . it is easy to obtain that $(2/1,2/2) = ^(2/1,2/2) = (2/1 “ exp( 0 . 000257 r)) 2 /i, 
/( 2 /i, 2 / 2 ) = (2/2, “ 2/1 +O.OOI2/2) and J(y) = (2/1, C(2/2))- 

Now, we will verify that system h3.1h\) defines aK— smooth discontinuous flow. First, condition (Cl) 
is verified since V<i>i( 2 /) = V$ 2 ( 2 /) = (1,0) 7 ^ 0, for all y G D. The jump function J{y) = (2/i,C(2/2)) 
continuously differentiable function. So, condition (C2) is valid. It is true that F fl F C (5F n i9F. In¬ 
equalities (V$i( 2 /),/( 2 /)) = ((1,0), ( 2 / 2 ,- 2/1 + 0 . 0012 / 2 )) = 2/2 7 ^ 0 and (V$ 2 ( 2 /),/( 2 /)) = ((1, 0), ( 2 / 2 ,- 2 /i + 
O.OOI 2 / 2 )) = 2/2 7 ^ 0, ify € F\9F, validate the condition (C4). Moreover, {V^i{y), f{y)) = ((1,0), (2/2,— 2 /i+ 
O.OOI 2 / 2 )) = 2/2 7 ^ 0 and (V$ 2 ( 2 /),/(//)) = ((1,0), ( 2 / 2 ,- 2 / 1 + 0 . 0012 / 2 )) = 2/2 7 ^ 0, if y Gr\df. Conditions 
(C 6 ) and (C7) hold as the function f is such defined. Thus, conditions (Cl) — (CIO) have been verified. 
Consequently, the system (|3.14l) defines the K— smooth discontinuous flow for all motions except the 
grazing ones. To complete the discussion, one need to linearize the system near the grazing solutions. 
First, we proceed with the linearization around the grazing periodic orbit iS.15\} . 

The solution, '!'(/) has two discontinuity moments 0 i = ^ and O 2 = a; in the interval [ 0 ,u;]. The 
corresponding discontinuity points are of ( 7 )— and (a)— types, respectively. Next, we will linearize the 
system at these points. The linearization at the second point exists w and the details of this will be 
analyzed in the next example. This time, we will focus on the grazing point y*. 

First, we assume that y(t) = y(t,0,y* + Ay), Ay = (Ai/i, A 2 / 2 ) is not a grazing solution. Moreover, 
the solution intersects the line Fi at time t = ^ near t = 9i as time increases. The meeting point 
y = (^ 1 ,^ 2 ) = ( 2/1 (C, 0, ( 2 /* + Ay)), 2 / 2 (/, 0 , ( 2 /* + Ay)), is transversal one. It is clear yi = exp(0.000257r) 
and y 2 > 0. In order to find a linearization at the moment t = 9i, we use formula (12.31) for y{t), and find 
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that 


dW,{y) 

dyl 


T(y) 


df{y{s))dy{s) ^^dT{y) ^ ^ ^ ^ j., , ^^dT{y) , ^ ^ ^ 

^ds + /(?/(s))^^ + Jyiyjiei + f[y{s))-^) + f{y{s) + J[y{s)))-;^ 


dy dy^ 


dy° 


9f{yis) + J{y{s))) dy{s) 


dx 


T-(v) 


dyi 


dVi 


ds, 


dy^ 

(3.17) 


where ei = (1,0)^, T denotes the transpose of a matrix. Substituting y = y to the formula (|3.17|) . we 
obtain that 


aw.mM +Ay)) ^ ^ ^,)) a^faK,o^g + A»)) 


dy° 


+Jy{y{^, 0, y* + Ay)) ( ei + f{y{f, 0, y* + j 




+/(y(e,0,(J(y(e,0,y* + A2/))))) 


dT{Jiyi^,0,y* + Ay))) 
dy° 


(3.18) 


Considering the formula (|3.10|) for the transversal point y = (j/i,?/ 2 ), the first component 


be evaluated as 


drjy) 


dT-jy) 

dy°i 


= — — . From the last equality, it is seen how the singularity appears at the grazing 
2/2 


point. Finally, we obtain that 


dW,{:y) 

dy<i 


2/2 

(-f) + 

1 0 

( ei + 

2/2 

-yi - O.OOlj /2 

^ 2/2 7 

0 -2Riy2 

V 

-yi - O.OOlj /2 


(3.19) 


-7?i(2/2)^ 

-yi + 0 . 001 Ri{y 2 )^ 



2/2 - 7?i(y2)^ 


1 0 


0 





yi + 0 . 0012/2 

-yi - 0 . 001(^2 - Riim)^) 

^ 2/2^ 

0 —2Riy2 




Calculating the righthand side of (13.1911 we have 


dWy{y) 

52/? 


- 7 ?i 2/2 - 1 

0.001(1 - i?iy 2 ) + 2i?i( 0 . 001^2 - yi) 


(3.20) 


The last expression demonstrates that the derivative is a continuous function of its arguments in a 
neighborhood of the grazing point. Since it is defined and continuous for the points, which are not from 
the grazing orbit by the last expression and for other points it can be determined by the limit procedure. 
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Indeed, one can easily show that the derivative at the grazing point y* is 


-1 


0.001- 1.8 exp(0.000257r) 


(3.21) 


Similarly, all other points of the grazing orbit can be discussed. 
Next, differentiating (ion with y[t) again we obtain that 


dW,{y) 

dy°2 


’■(y) 




dy dy'i 


dvl 


dy'i 


df{y{s) + J{y{s))) dy{s) 


dx 


’-(v) 


dy° 


ds. 


(3.22) 


where e^ = (0,1)^. Calculate the right hand side of (13.221) at the point y = (yi, j/ 2 ) to obtain 


3W.iyi(M + Ay)) ^ ^ ^^^^ ar(y((.0 y- + Ay)) 


dy° 


NJyiyi.^, 0 , y* + Ay)) 62 + /(y(c, 0, y* + Ay))) 


To calculate the fraction 


+/(y(C,0,y* +Ay)) 
5r(y(C,0,y* + Ay)) 


dy^ 

dT{y{^,0,y* + Ay))^ 


dy^ 


dT{y{^,0,y* + Ay)) 

ayO 


(3.23) 


5y2° 


in (I3.23p . we apply formula (13.101) for the transversal 


^t(v) Ot(v) 

point y = (yi,y 2 ). The second component ^ takes the form ^ = 0. This and formula (13.231) 

oy2 oy2 

imply 


aWi(y) 

5y2° 


0 

-2i?y2 


(3.24) 


Similar to (|3.21|) . one can obtain that 


dW,{y*) 

dy°2 


0 

0 


(3.25) 


Joining (13.211) and (13.251) . it can be obtained that 


W.y{y*) 


-1 0 
0.001-1.8 exp(0.000257r) 0 


(3.26) 
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The continuity of the derivatives in a neighborhood of y* implies that the function W is differentiable at 
the grazing point y = y*, and the condition (A2) is valid. 

Now, on the basis of the discussion made above, one can obtain the bivalued matrix of coefficients for 
the grazing point as 


O 


2 , 



-1 


0.001 - 1.8exp(0.000257r) 


0 

0 


if (Nl) is valid, 
if (N2) is valid. 


The matrix = O 2 is for near solutions of ic).15\} which are in the region where Z in, see Fig. QJ 
and do not intersect the curve of discontinuity T 1 . The matrix 



-1 

0.001- 1.8 exp(0.000257r) 


0 

0 


is for near solutions of iS.15\} . which intersects the curve of discontinuity Fi. They start in the subregion, 
where the point Y is placed. Thus, the linearization for 'I' (i) at the grazing point exists. Moreover, since 
another point of discontinuity (0, exp(0.00057r)) is not grazing, the linearization at the point exist as well 
as linearization at points of continuity Conseguently, there exist linearization around 'l'(t). 

To verify condition (A3), consider a near solution y(t) = y(t,0,y) to '^{t), where y = ( 0 ,^ 2 ), y 2 > 
'I'2(0) = 1, which satisfy the condition (iVl). It is true that ~ “ '2 “ "2 ’ first coordinate of 

the near solution is yi(t) = y exp(0.0005i) sin(t) and 

yi{^) = yi('^) = y exp(0.000257r) > exp(0.000257r) = —). Thus, the meeting moment of near solution 

yit) with the surface of discontinuity is less than So, it implies that 0 < r(y) < § ~ e for a small 
number e if the first coordinate ofy is close to exp(0.000257r). This validates condition (A3). Now, Lemma 
ro proves the condition (C). 

Now, let us consider the point (0, 0). We have that (V<I>((0,0)), /((0,0))) = ((1, 0), (0, 0)) = 0. That is 
the origin is a grazing point. In the same time it is a fixed point of the system. For this particular grazing 
point, we can find the linearization directly. Indeed, all the near solutions satisfy the linear impulsive 
system, 

x'l = X2, 

x '2 = —xi + 0.00112, (3.27) 

Aa;2Ui=o = -(1 + R2)x2- 

Consider a solution x(t) = x{t,0,xo), where Xq = (x^jX^) 7 ^ (0,0) with moments of discontinuity 
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9i^i € Z, then the linearization system for the equation around the equilibrium is 


Ml = M2, 


M2 = —Ml + O.OOIU 2 , 


(3.28) 


^U2\t=Si = -(1 + R2)u2. 


Indeed, ifui{t), mi(0) = ei, U 2 {t), M2(0) = 62, are solutions of (13.281) . then one can see that xft) — (0,0) = 
XiUi{t) +X 2 U 2 it), for all t G K. 

We have obtained that linearization exists for both grazing solutions '^(t), and the equilibrium at 
the origin. Moreover, conditions (Cl) — (CIO) are valid and all other solutions are B-differentiable in 
parameters m- Thus, the system (13.141) defines a K— smooth discontinuous flow in the plane. 

In the next example, we will finalize the linearization around the grazing solution 4'(t). 

Example 3.2 (Linearization around the grazing discontinuous cycle). We continue analysis of the last 
example, and complete the variational system for 'l'(t). 

Let us consider this time, the linearization at the non-grazing moment w = tt. The discontinuity point 
is c = (0, — exp(0.00057r)) and it is of (a)— type, since 


(V$(c),/(c)) = ((l,0)(-exp(0.00057r),-0.001 exp(0.00057r))) = - exp(0.00057r) 7 ^ 0. 


By using formula itlKA) . one can compute the gradient as Vt^c) = (exp(—O.OOOS tt), 0). 


Then, utilizing S/rfc) and formula L‘L11\) . one can determine that the matrix of linearization at the 


moment tt is 


exp(—O.OOOS tt) 0 


0.001 


0 


From the monotonicity of the jump function, —Riy^, it follows that the the yellow and blue subregions 
of G are invariant. Consequently, for each solution near to 'l'(t), the sequences Bi is of two types 


Bi = i G Z and j = 1,2, where i^ 2 i-i = R^ 2 i-i — 


-1 


0 


0.001-1.8 exp(0.000257r) 0 



, * G Z. That is, the condition (A4) is valid and the linearization 


0.001 0 


around the periodic solution i3.15\} on R is of two subsystems: 
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Ml = U2 


u'2 = —Ml + O.OOIM2, 

AM|t=e2i = 


(3.29) 


and 


Ml = M2, 

M2 = —Ml + O.OOIM2, 

(3.30) 

AM|t=e2._i = 

AM|t=e2, = 

where 02 i-i = and 02 i = 

r/ie sequences {Dp'*}, } = 1,2, are 2— periodic. It is appearant that system iS.29\} + iS.SfXl is a 
(w, 2)— periodic. Thus, the variational system for the grazing solution is constructed. 


4 Orbital stability 

In this section, we proceed investigation of the grazing periodic solution 'l'(t). Analysis of orbital stability 
will be taken into account. Denote by B{z,d), an open ball with center at z and the radius (5 > 0 for 
a fixed point z G T \ 9r. By condition (C3), the ball is divided by surface F into two connected open 
regions. Denote c“''(z, (5), for the region, where solution x{t) = x{t, 0, z) of (12.21) enters as time increases. 
The region is depicted in Figure [31 

Set the path of the periodic solution 'l'(t) as 


p ■.= {x € D : X = 'l'(t), t G K}. 



Define dist{A, a) = infag^i ||q: — a||, where A is a set, and a is a point. 
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Definition 4.1 The periodic solution 'f'(f) : R —>■ D o/ i2.1\) is said to be orbitally stable if for every e > 
0, there corresponds 6 = (5(e) > 0 such that dist{x{t^0,xo)^r]) < e, for all t>0, provided dist(xo,r]) < S 
and xq ^ UiC+(4'(0i), (5), for * = 1,... ,m, where m is the number of points 4'(0i) G r\ 9r. 

The point Xq is not considered in regions 5), i = 1,..., m, since solutions which start there 

move continuously on a finite interval, while 'h (t) experiences a non-zero jump at t = Oi and this violates 
the continuity in initial value, in general. In the same time, we take into account any region adjoint to 
points of (9r, since the jump of d'(f) is zero there and, consequently, the continuous dependence in initial 
value is valid for all near points. 

Definition 4.2 The solution 'l'(t) : D of i2.1\) is said to have asymptotic phase property if a 

(5 > 0 exists such that to each Xg satisfying dist(xg^r]) < 5 and xq ^ UiC“'"('I'(0i),5), for i = l,...,m, 
there corresponds an asymptotic phase a{xo) € R with property: for all e > 0, there exists T{e) > 0, such 
that -f Qf(a;o), 0, a:o) is in e-neighborhood of^(t) in B—topology for t £ [T(e),oo). 

Let us consider the following system, which will be needed in the following lemmas and theorem 


x' = A{t)x, 
Ax|t=^, = BiU, 


(4.31) 


where A{t) and Bi are nx n function-matrices, A{t-\-ui) = A{t), for alH € R and there exists an integer 
p such that = Q oj and Bi^p = Bi, for all 5 G Z. 


Lemma 4.1 Assume that system {4-31^ has a simple unit characteristic multiplier and the remaining 
n—\ ones are in modulus less than unity. Then, the system {4-^^^ has a real fundamental matrix X{t), 
of the form 


X{t) = P{t) 


(4.32) 


V 


0 


exp [Bt) 


where P £ PC'^(R, 0) is a regular, oj-periodic matrix, and B is an (n — 1) x (n — 1) matrix with all 
eigenvalues have negative real parts. 


Proof. Denote the matrix A(t),A(0) = I, as fundamental matrix of system (14.311) . There exists a 
matrix Bi such that the substitution x = P{t)z, where P{t) = X{t) exp(—Pit), transforms (I4.31|) to the 
following system with constant coefficient [T], 


z' = Az. 


(4.33) 


The matrix exp(Aa;) has a simple unit eigenvalue and remaining (n — 1) ones are in modulus less 
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than unity. Hence, there exists real nonsingular matrix M, which satisfies 


M ^ exp(Aaj)M 


1 0 
0 Cl 


The remaining part of the proof is same as proof of Lemma 5.1.1 in HZ]. □ 

Throughout this section, we will assume that (H4) is valid. That is, the variational system (13.131) 
consists of m periodic subsystems. For each of these systems, we find the matrix of monodromy, Ujioj) 
and denote corresponding Floquet multipliers by i = 1,..., n, j = 1,..., m. In the next part of the 
paper, the following assumption is needed. 

(A5) = 1 and | < 1, i = 2,..., n for each j = 1,..., m. 

Lemma 4.2 Assume that the assumptions (A4) and (A5) are valid. Then, for each j = 1,... ,m, the 
system i&.lA) admits a fundamental matrix of the form 


Uj{t) = Pj{t)[l,exp{HjUj)], t G R, (4.34) 

where Pj G PC^(R, C) is a regular, w-periodic matrix and Hj is an (n — 1) x (n — 1)— matrix with all 
eigenvalues have negative real parts. 

The proof of Lemma 14.21 can be done similar to that of Lemma 14.11 

Theorem 4.1 Assume that conditions (Cl) —(C7), (CIO), and the assumptions (Al) —(A5) hold. Then 
UJ- periodic solution 4'(t) of (12.11) is orbitally asymptotically stable and has the asymptotic phase property. 

Proof. Since of the group property, we may assume 'I'(O) is not a discontinuity point. Then, one can 
displace the origin to the point 4'(0), and the coordinate system can be rotated in such a way that the 
tangent vector 'I'q = 4''(0) points in the direction of the positive x\ axis i.e. the coordinates of this 
vector are 'I'q = ('I'q^, 0,..., 0), > 0. 

Let 0i, i G Z, be the discontinuity moments of 'l'(t). Denote the path of the solution hy rj = {x G X : 
X = 'it{t),t G K}. There exists a natural number p, such that 9i+p = 9i + ui for all i. Because of condi¬ 
tions (Cl) — (C7) and AT—differentiability of 'l'(t) there exists continuous dependence on initial data and 
consequently there exists a neighborhood of rj such that any solutions which starts in the set will have 
moments of discontinuity which constitute a B— sequence with difference between neighbors approxi¬ 
mately equal to the distance between corresponding neighbor moments of discontinuity of the periodic 
solution Consequently we can determine variational system for 'l'(t), with points of discontinuity 

di, i G Z. 

On the basis of discussion in Section 2.1, one can define in the neighborhood of p a B— equivalent 
system of type (1^ . The variational system of it takes the form 
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z = A{t)z + r{t, z), 

=n^^^z + qi(z), j = l, 2 ,...,m, 


( 4 . 35 ) 


wheier(t,z') = [f +z) - - A(t)z and q,(z) = Wi(d>(Oi) + z)-W.id'id,)) - z, are contin¬ 
uous functions, and matrices satisfy condition {AA). The functions are continuously differentiable 

with respect to 2 ;. One can verify that r(t, 0) = qi{A) = 0 and r{t + u!,z) = r{t, z) for t S K.. Moreover, 
the derivatives satisfy r'{t,0) = ^'^(O) = 0 and the functions r{t,z) —>• 0 , qi{z) —>■ 0 , r'^{t,z) —>■ 0 and 
q'i^iz) —>■ 0, as z 0 uniformly in t € [0, 00 ), i> 0. Each system (14.351) for j = 1,2,... ,m, corresponds 
to a region adjoint to initial value, xq such that these regions cover a neighborhood of xq. 

Fix a number j and denote Yj(t) the fundamental matrix of adjoint to (14.351) linear homogeneous 
system 


y = Mt)y^ 

Ay\t=e^ = 


( 4 . 36 ) 


of the form ( 14 . 341 ) . One can verify that 


Yjit)Y-\s) = P,{t) 


^ 1 0 

^ 0 exp(iJj(t — s)) 




for —00 < t, s < 00 . 
We can write 


( 4 . 37 ) 


/ 


/ 


^ 0 exp{Hj{t — s)) J y 0 exp(i 7 j(t — s)) 
where 0 „_i is the (n — 1 ) x (n — 1 ) zero matrix. Then it can be driven 


1 0 
0 0„_i 


Yj{t)Y-\s) = G[Y{t, s) + s) = G^Yit, s), 

where 


G«(t,s) = P,(t) 


0 0 
0 exp{Hj{t — s)) 

/ 




Gi^\t,s)=P,{t) 


1 0 
0 o„_i 


P-\S). 
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Denote the eigenvalues of the matrix Hj by ■ ■ ■, Xn^ ■ By means of the Lemma HTT] and 1?^ there 
exits a number a > 0, such that Re{X^^^) < —a, fc = 2,3,..., n, where Re{z) means the real part of the 
number, z. Taking into account that the matrices Pj and P~^ are regular and periodic, the following 
estimates can be calculated 

s)| < exp{—a{t — s)), (4.38) 

\Gi^\t,s)\<K^^\ (4.39) 

where is a positive real constant. 

Denote the first column of the fundamental matrix Y by By the equation (I4.34|) . is equal to 
the first column of Pj , this means that it is a w-periodic solution of (13.131) . 

By assumptions of the theorem the variational system (I4.35|l satisfies the conditions of Lemma 14.21 
and one can verify that the following estimate is true m 

\Yj{t)\ < k[^'' exp{—at) for t>0, (4.40) 

where is a positive constant. Let us setup the following integral equation 

t OO 

{t,a) = Yj{t)a + j G^^\t,s)r{s,z{s))ds — J G 2 \t, s)r{s, z{s))ds 

0 t 

+ ^ G[^\t,9k+)qk{z{0k)) — ^ G2\t,9k+)qk{z{dk)), (4-41) 

0 <dk<.t t<dk<oo 

where a = [0, 02 ,..., a„], G M, i = 2,3, ...,n, are orthogonal to d''(0), i.e. with the zero first 
coordinate. 

Let z^\t,a) = 0, and consider the following successive approximations 

00 

^^■^^(^, 0 ) = Y,(t)a+ / G^^\t,s)r{s,Zk-i{s))ds+ ^G^^\t,9k+)qk{zk-i{dk)), (4.42) 

0 ^==1 

for k = 1,2,.... By using the approximation (14.421) and estimation (14.401) . one can verify that 

\z[^\t,a)\ < K[^^\a\exp{—at/2). (4.43) 

We will show that the bounded solution of (14.411) exists and satisfies (14.351) . For arbitrary positive small 
number L, there exists a number S = S{L) such that for \zi\ < 5, \z 2 \ < d 

|r(f, zi) - r{t, Z 2 )\ < L\zi - Z 2 l (4.44) 
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and 


\qi{zi) - qi{z 2 )\ < L\zi - Z 2 \, (4.45) 

uniformly in i € [0,c»). 

Denote by Li = (— — - -^• 

Va l-exp{-ae/2)J 

Next, by using mathematical induction, we are going to show that zi^\t, a), s = 1, 2,..., are defined 
for t £ [0, oo) and satisfy 

kl+i(i,a) - < itrp Vl exp(-at/2)/2®, s = 0,1,2,..., (4.46) 

if L < Li. Utilizing Lemma [4.21 and inequalities ()4.40|1 . (14.4411 . (14.451) and 9i+i — 9i > 9_,i £ Z, one can 
verify that 

\zi^lT^{t,a) - zi^\t,a)\ < K[^'^\a\Liexp{-at/2)/{2'^a). (4.47) 

As a consequence of (I4.46L the sequence z^'^^(t,a) converges uniformly on t £ [0,oo), |a| < S/2K[^\ 
and 

\zi^\t,a)\ < 2K[^'>\a\exp{-at/2),s = 1,2,.... 

Therefore, the limit function exists on the same domain, it is piecewise continuous, satisfies 

(14.411) and the following estimate 


|z^'^^(t,a)| < 2Arp^|a| exp(—Q;t/2). (4.48) 

Denote by z{t) = z^^'> {t, a), for j = 1,2,... ,m. Next, we will verify that z^^'^ {t, a) satisfies (14.351) . For 
it, differentiate ()4.41|) 

t 

z\t) =Y^{t)a^G^^\t,t)r{t,z{t))+G{(t,t)r{t,z{t))+ J G[\\t,s)r{s, z{s))ds 

0 

oo 

- f Gil\t,s)r{s,z{s))ds + ^ G[\\t,9k+)qk{z{9k)) - ^ G^2iHt^^k+)qk{z{9k)) 

^ 0 < 9 k<t t< 9 k<oo 

oo 

= A{t)Yj{t)a + G^^'>{t,t)r{t,z{t)) + f A{t)G^^'> {t, s)r{s, z{s))ds + ^ A{t)G^^'> {t,9k+)qk{z{9k)) 

g 0<ei<t 

= A{t)z{t) + r{t, z{t)). 

Fix 0fc, /c £ Z, then 
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z{Ok+) - z{9k) 


dk oo 

= Yj {Ok+)a + J G[^'’ {0k+, s)r{s, z{s))ds - J (0^+, s)r(s, z{s))ds 

0 Sfc 

+ G^i^ {9k+,0i+)qi{z{9i+)) — G^2^ (9k+,9i+)qi(z{9i+)) 

0<9i<9k 6k<0i<oQ 

9k oo 

-Yj{9k)a- J G[^\9k,s)r{s,z{s))ds + J G^^^dk, s)r{s, z{s))ds 

0 Sk 

— G^i^ {9k,9i+)qi{z{9i)) + G2^ {9kT9i+)qi{z{9i+)) 

O<0x<0fc 0fc<0x<OO 

= D''^'>z[9k)+qk{z{9k)). 

The above discussion proves that {t,a),j = 1,, 2,..., m, are bounded solutions of system (14.351) . 

We will determine the initial values of bounded solutions in terms of (n— 1) parameters , ■ • ■, a™ ^, j = 
1,, 2,..., m. Denote = [q, a^, Og,..., a^]. By using (14.411) . we obtain 


=Yj{0)a^^'^ - f G^2\0,s)r{s,z{s))ds - ^ G^2\o,9k+)qk{z{9k)) 

n 0<@fc<oo 



In the way utilized in one can show that the coordinates of the initial value {xi ,..., x„) £ D of the 
solution z^Y satisfy the equation 

n 

Y'^c’iXi - hj{x 2 , ...,Xn)=0, (4.49) 

i^2 

where hj £ C^,j = 1,, 2,..., m. 

One can see that equation (14.491) determines (n—1) dimensional hypersurfaces C D, j = 1,2,... ,m, 
in a neighborhood of the origin such that each solution which starts at the surface satisfies inequality 
(|4.48l) . From the analytical representation, it follows that the equation of the tangent space of at the 

n 

origin is described by the equation xi + ^ and the first coordinate of the gradient of the left hand 

i=2 

side in (I4.49|) is unity. Moreover, the path -q intersects transversely. This and condition {AA) imply 
that the path of every solution near 4'(t) intersects one of the manifolds SGj = 1, 2,..., m, at some 
t £ [0, 2ujj . 

Because of the continuous dependence on initial values, a S{e) > 0 exists for a given e > 0, such that 
if dist{x^,qs) < ^{e), then the solution (j){t,x^) is defined on [0,2a;], and dist{(j){t,x^),q) < e < ei for 
t £ [0, 2T]. Therefore, the path of (j){t,x^) intersects for some j = l,2,...,m and ti £ [0,2a;]. The 
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solution ,x^)) = (f>{t + has its initial value in , consequently, satisfies (I4.48|) . In the light 

of the B— equivalence, the corresponding solution x(t), a:(0) = 0(0) —4'(0), of (14.351) satisfies the property 
that for all e > 0, there exists T(e) such that x{t) is in an e— neighborhood of 4>(t) for t G [r(e),oo). 
That is, the solution \E'(t) is orbitally asymptotically stable and there exists an asymptotical phase. □ 

Definitions of the orbital stability and an asymptotic phase as well as theorem of orbital stability for 
non-grazing periodic solutions are also presented in |44| . In our paper, we suggest the orbital stability 
theorem for grazing periodic solutions, its proof and formulate the definitions for the stability. They are 
different in many aspects from those provided in [44]. It is valuable that they also valid, if the solution 
is non-grazing. 

To shed light on our theoretical results, we will present the following examples. 

Example 4.1 We continue with the system presented in Examvles \3.1\ and \3.‘A In Examvle \3.1[ we 
verified that system iS.16\} defines a K— smooth discontinuous flow in the plane and the variational 
system h3.2fA) + h3.8(Kl around the grazing periodic solution, 'l'(t) is approved. 

Using systems i3.291) and i3.30\} . one can evaluate the Floquet multipliers as = 1, = 0.8551, 

= 1 and = 0. This verifies condition (A5). 

The conditions ((71) — ((77) and ((710) are validated and the assumptions (A4-) and (A5) verified. By 
using Theorem E3 we can assert that the solution, 'l'(t) is orbitally asymptotically stable. The stability 
is illustrated in Fig. The red one is for a trajectory of the discontinuous periodic solution i3.15\) of 
i3.14\ l and the blue one is for the near solution of with initial value yo = (0.8,1.2). It can be 

observed from Fig. that the blue trajectory approaches the red one as time increases. 



Figure 4: The red discontinuous cycle of (13.1411 axially grazes F at (0.000257r), 0) and (0, — exp(—0.00057r)) 
is an (Q:)-type point. The blue arcs are of the trajectory with initial value (0.8,1.2). It can be observed 
that it approaches the grazing one as time increases. 


Example 4.2 (A periodic solution with a non-axial grazing). We will take into account the following 
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autonomous system with variable moments of impulses 


x'l = X2, 


x'2 = -xi, 

Ax2|a:er = 


1 

72 

1 


Xi + K{x 2 — Xif, 


X2 + K{x 2 - Xif, 


(4.50) 


where F = {(xi^X 2 )\xi + X 2 = V^}, f = {{xi^X 2 )\xi = X 2 } and K = 0.11. It is easy to verify that the 
point X* = ( 2 ’ 2 ^ ^ grazing point because {\7 t[x*), f{x*)) = ((1,1), (i, —i)) = 0 and the grazing is 

non-axial. We assume that the domain is the plane. 

The solution '^(t) = (sin(i), cos(i)), t G M. is a grazing one, since the point x* = is from its 

orbit. The cycle and the line of discontinuity are depicted in Figure O 



Figure 5: The red curve is the orbit of 'l'(t) which grazes non-axially the line of discontinuity. 

Let us consider the linearization at the grazing point x* next. We will consider the near solution 
x(t) = x(t,0,x* + Ax). Denote t = the moment when the solution meets the surface of discontinuity 
r at the point x = x{ff) = x(^, 0, x* + Acc). Taking into account formulae (I3.17|l . (111.181) with (I4.5()|) . one 
can obtain the following matrix 


dWi{x{£,, 0, X* + Ax)) 


X2 

1 

_ + 

— 2K{x2—Xi) —2K{x2 — Xi) 

-Xi 

Xi — X2 

— 2K{x2—Xi) —2K{x2 — Xi) 


X ei + 


X2 


-^+K{x2-Xi)^ 

-Xl 

XI-X 2 I 

-K{x2-Xi)^ 


Xi — X2 


(4.51) 


Calculating the right hand side of the expression (14.5111 . we obtain that 
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(4.52) 


dWi{x{^,0,x* + Ax)) 


-y/2 + 0.22' 

2V^ 

^2 + 0.22 
2^2 


Using similar method with that of the first one, the second derivative can be computed as 


dWi{x{f, 0, X* + Ax)) 

dx^ 


V2 + O.22 
2 V 2 

-^2 + 0.22 
2V^ . 


(4.53) 


Combining (14.521) and (14.531) . we can obtain the following matrix for the linearization at the grazing 
point X *, 


W,,,{x*) = 


- y /2 + 0.22 ^2 + 0.22 

2^2 2V2 

^2 + 0.22 -V2 + O.22 


2^2 


2y/2 


(4.54) 


It is appearant that the matrix Wix{x*) is continuous with respect to its arguments, since it is constant if 
the point x* + Ax is not from the orbit of the grazing solution. Since of the limit procedure, it is the same 
constant for all points of the grazing solution. Thus, the Jacobian is constant matrix in a neighborhood 
of the grazing point and condition {A2) is valid. 

Now, let us check the validity of the condition (A3). Consider a near solution x(t) = x(<:,0,x), to 
the grazing cycle '^{t), where x = (0,X2), X 2 > '£' 2 ( 0 ) = 1. So, the near solution x{f) satisfies the 

condition {Nl). For the grazing periodic solution, it is true that 9i+i —9i =2 'k = oj. The grazing solution 

'£' [t) = x(t, 0, (0,1)), touches the line of discontinuity T att=^. The first coordinate of the near solution 
is Xi{t) = X 2 sin(t), and Xi( —) = X 2 sin( —) = —^ > '£'i( —) = —p. Consequently, the near solution x{t) 

8 8 ^J2 8 ^J2 

meets the line of discontinuity T before the moment j. This implies that 0 < t(x) < j — e, for a 
small positive e whenever xi{t) is close to ■^. Thus, the condition (A3) is valid and Lemma \S.l\ vroves 
condition (C). 

In the light of the above discussion, the bivalued matrix of coefficients for the grazing point is easily 
obtained as 
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if (N1) is valid, 


O 


2 , 


[-^2 + 0.22 72 + 0.22' 
272 272 

72 + 0.22 -72 + 0.22 

. L 275 27I - 


if (N2) is valid. 


(4.55) 


It is appearant that the interior of the grazing orbit is invariant. Let us show that the external part 
of the unit circle is positively invariant. It is sufficient to demonstrate that Ji(a:i)^ + J 2 {x 2 )'^ > 1 for 
any {xx,X 2 ) S F. Denote xi = z and X 2 = 72 — 2 ; and consider the formula 


^ + 0.11 

72 


where F{-^) = 1- It is easy to calculate that = 0 o,nd F"{^) = 0.8872 > 0. Consequently, 


Fi.) - F(^) = f (.) - 1 = - +)^ + „(||. - +f) > 0. 


if z is close to . Thus, near the grazing point, the external region is invariant. From this discussion, 
since of the formula (I4.55|l . we can conclude that the condition {AA) is valid. Taking into account it with 
the expression (I4.55|) . the linearization system for (|4.50|) around the grazing solution 'F(i) is obtained as 


(72 - 2z)Y + (^ + 0 . 11(72 - 2z)y 
72 


F{z) = Mzf + J2{\/2 - zf = ( 


u'l = U2, 

u '2 = —U\, 

Au(27ri) = d \^^ u , 


(4.56) 


where = O 2 and = 


-72 + 0.22 72 + 0 . 22 ' 

272 272 

72 + 0.22 -72 + 0.22 


i e z. 


272 272 

To finalize stability analysis, consider the first system in (14.561) . with matrices = O 2 . Its multi¬ 
pliers are = 1 and it constitutes the linearization for the orbits which are inside the circle. The 

system does not give a decision by orbital stability theorem. Theorem EH Nevertheless, from the simple 

analysis m result, we know that the grazing orbit is stable with respect to inside orbits of the system. 

( 2 ,') ( 2 ) 

The linearization of orbits which are outside of the circle has multipliers p7 = 1 and P 2 = —0.15. It 
means that the periodic solution is orbitally stable with respect to solutions outside of the circle. Summa¬ 
rizing the discussion, we can conclude that the periodic solution is stable. The stability result is observed 
through simulations and it is seen in Fig. 0 
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^ 0 
- 0.5 

-1 

- 1 . 5 '--‘‘-■-‘- 1 

- 1.5 -1 - 0.5 0 0.5 1 1.5 

Xi(t) 

Figure 6: The red orbit of system (|4.50p non-axially grazes the surface T. The magenta trajectory with 
initial point (0,1.32) approaches the cycle as time increases. The green cycle with initial point (0, 0.96) 
demonstrates the inside stability of the grazing orbit. 

5 Small parameter analysis and grazing bifurcation 

In this part, we will discuss existence and bifurcation of cycles for perturbed systems, if the generating 
one admits a grazing periodic solution. In continuous dynamical systems, a small parameter may cause a 
change in the number of periodic solutions in critical cases. In the present analysis, we will demonstrate 
that the change may happen in non-critical cases, since of the non-transversality. That is why, one can 
say that grazing bifurcation is under discussion. Let us deal with the following system 



x' = fix) + figix,g.), 
Aa;Ugr(A.) = /i). 


(5.57) 


where x € € M, r(p) = {a;| ^{x)+fj,<j){x, fi) = 0}, /r G (—Mo, Mo), and mo is a sufficiently small positive 

number. Functions/(x),/(x) and <I>(a:) are continuously differentiable up to second order, ^(a;, m), 7L(a;, m) 
are continuously differentiable in x and fi. The function (fix, g) is continuously differentiable in x up to 
second order and to first order in m- We assume that the generating system for (15.571) is the system (12.11) 
with all conditions assumed for the system, earlier. The main assumption of this section is that dUD 
admits a w—periodic solution, 'l'(t). Let 'I'(O) = (C°, C 2 , • ■ •, Cn) be the initial value of the solution. 

Our aim is to find conditions that verify the existence of periodic solutions of (15.571) with a period 
such that for m = 0, the periodic solutions of (15.571) are turned down to 'l'(t). It is common for the 
autonomous systems that the period does not coincide with uj. Thus, in the remaining part of the 
paper, we will consider the period as an unknown variable. 

Since ^'(0) is not an equilibrium, there is a number j = 1, 2,.. ., n, such that /j(Ci, Ci, ■ • ■, Cn) 7^ Cl¬ 
in other words, the vector field is transversal to line Xj = near the point. Hence, to try points near to 
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'I'(O) for the periodicity, it is sufficient to consider those with j—th coordinate is equal to [31]. For the 
discontinuous dynamics, the choice of the fixed coordinate can be made easier if the surface of disconti¬ 
nuity is provided with a constant coordinate. We will demonstrate this in examples. Denote the initial 
values of the intended periodic solution by ^i, C 2 ,..., Cn- Assume that one initial value Q is known, i.e. ^1*. 
Thus, the problem contains n—many unknowns, they can be presented as Ci, C 2 , ..., Ci-ij Cj-i-ij • ■ • i Cm 
Denote the solution of (|5.57l) by Xs(t, (^ 1 X 2 , ■ ■ ■, Cn, p) with initial conditions a;s(0, Ci, ^ 2 ,..., Cn, p) = Cs- 
To determine the unknowns, we will consider the Poincare criterion, which can be written as 


^fe(^,Cl,C2,- ■ ■ ,Cn,p) = Xki.^,ClX2,---Xn,X - Cfc = 0, fc = 1, 2, . . . , n. 


(5.58) 


where Q = (j. The equations (15.5811 are satisfied with fi = 0, ^ = w, Ci = C?, i = 1, 2,..., n, since 
is the periodic solution. 

The following condition for the determinant is also needed in the remaining part paper. 

(A6) 



.0 .0 



.0 .0 

-.C^o)) 




a{,y‘2{^,ClX°2^ 

aCn 



asr 



aCn 



1 


a(j^„(ojx?Xl 




asr 



aCn 



0 (5.59) 


Theorem 5.1 Assume that condition (A6) is valid. Then, 15.571) admits a non-trivial periodic solution, 
which converges in the B— topology to the non-trivial lo- periodic solution of 15.571) as /i tends to zero. 

We will present the following examples to realize our theoretical results. 

Example 5.1 In this example, we will consider the perturbed system in case the generating system has 
a graziness. To show that, let us take into account the following perturbed system 


(5.60) 


Xi = X2, 

x '2 = — 0.00112 — Xi, 

Aa;2UGri = -(1 + R 1 X 2 + fJ-X2)x2, 

Aa; 2 |xGr 2 = -(1 + R 2 + Kx 2 - exp(0.00l7r/2))x2. 

It is easy to see that the system 15.6'dl) is of the form /15.57I) . For p = 0, the generating system became 
For the perturbed system /15.6‘dl) . we will investigate existence of the periodic solution around the 
grazing periodic solution of with the help of Theorem 15.11 

There are two sorts of possible periodic solutions of 15. 60\} around the grazing one. One of them 
has two impulse moments during the period since it crosses both lines of discontinuity, i.e. Xi = 0 
and xi = exp(0.000257r). The other sort is the periodic solution which does not intersect the line Xi = 
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exp(0.000257r) and intersects the line xi = 0. We will show the existence of both type of periodic solutions 
if \fi\ sufficiently small. 

Let us start with the second type, assume that the solution for the perturbed system exists and it starts 
at the point (0, a;o 2 ), a:o 2 < 1 and does not intersect the line xi = exp(0.000257r). Denote the initial values 
of the periodic solution by and (^ 2 - Since the periodic solution necessarily intersects the line xi = 0, 
one can choose Ci = Cj* = 0. By specifying the formula in (|5.58D for the system (I5.60p . it is easy to obtain 
the following expressions 


^i(^,0,C2,/i)=a;i(=;^,0,C2,M)=0, 

0, C2, m) = X2(^, 0, C2, m) - C2 = 0. 


(5.61) 


Next, taking the derivative of the expressions in (15.611) . we can obtain the following 


a{.5^i{3r,0X2w)) d{j^i{3r,0X2,v)) 

dsr dC2 

diJ^2{S,0X2w)) d{5^2{3r ftX2,u)) 

dsr dC2 

The determinant ()5.62l) is calculated by means of the monodromy matrix of with the impulse 

matrix = O 2 , i-e. 


dxi (tj.O.Cs.O) 

dS 

3x2((.iJ,0,C2 , 0 ) 

as 


0a:i(u;,0,C2 .0) 

aC2 

ax2{i.ii,oxfo) 

a(2 


- 1 


(5.62) 


1 -0.0317 

1.0158 -0.1014 


(5.63) 


Taking into account the system (15.661) with (15.631) at C 2 
that 


(2 and SL = u} for ^ = 0, one can derive 


as 

aS’2{u2XX2.X) 

as 


0 C 2 

as2(oj,oxfo) 1 

a<2 ^ 


-0.0317 exp(0.000257r) f 0. 


(5.64) 


This verifies condition (^46). Thus, condition (^46) is valid, then by utilizing Theorem (15.11) . we can 
assert that the system if 5. 571) admits a non-trivial periodic solution, which converges in the B— topology 
to the non-trivial uj-periodic solution of a as pL tends to zero. 

Now, let us verify that system h5.6(A) has a circle which intersects the line xi = exp(0.000257r) in the 
neighborhood 0 /(exp(0.000257r), 0). So, the periodic solution will attain two discontinuity moments in a 
period. Denote the initial values of the periodic solution by and C 2 . To apply the condition {A6), fix 
one initial value Ci = Ci = 0 0 / the intended periodic solution and in the light of the expressions (15.581) 
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^i(^,0,C2,/i)=a;i(=;^,0,C2,M)=0, 

^2(^, 0 , C2, m) = X2(^, 0 , C2, m) - C2 = 0 . 


(5.65) 


Taking the derivative of the expressions (15.651) with respect to variables ^ and ^2 j one can obtain the 
following 


a(.5"2(5",0.C2.M)) 

a(j5^i(j7,o,C2.M)) 


dxi (a;,0,Cj,0) 

dxi{ujXX 9 ., 0 ) 

dsr 

SC 2 


asr 

ac2 

a(j>^2(J7.o.C2.M)) 

d(y2{3’,^,Q2 ,u)) 


dx 2 {uj, 0 X 9 , 0 ) 

dx2(i.Li,0,C2,0) 1 

dsr 

dC2 


as 

a(2 


To determine the above determinant, the monodromy matrix of i3.1S\) with the jump matrix can 
be evaluated as 


1 0.01 

0 0.704 

For p, = 0, with the values w and (2 the determinant (15.661) can be determined as 


(5.67) 


0(jy2(<.^,o,cg.o)) 

dsr 


a(.yi(c^,o,Cg.o)) 

0 C 2 

0 X 2 , 0 )) 

9C2 


0 0.01 
- exp(0.000257r) -0.296 


0.01 exp(0.000257r) ^ 0. (5.68) 


This verifies condition (A6). So, By Theorem \5.1[ we can conclude that the perturbed system (15.601) 
admits a non-trivial periodic solution which converges in the B— topology to the non-trivial 

uj-periodic solution of jS.lfl ) as p, tends to zero such that ^(0) = uj. 

In Fig. [3 some numerical results are provided to show the solutions of system S5.60i} with p = 0.05. 

The periodic solutions for /i ^ 0 are not grazing. For fi = 0, we have one periodic solution which 
is orbitally stable, and for fj, < 0, there exist two periodic solutions. One of them has one discontinuity 
moment in each period, in other words, the cycle does not intersect the surface of discontinuity around 
grazing point and it is orbitally stable and the other one has two discontinuity moments in each period. 
This means, the number of periodic solutions increases by variation of pL, around fi = 0. So, we will call 
that bifurcation of periodic solution from a grazing cycle. 


Example 5.2 Let us consider the following system with variable moments of impulses and a small 
parameter 

x\ = X2, 

x '2 = -0.0001[x2 + (cci - 1)^ - (1 + h)'^]x 2 - a:i + 1, (5.69) 

^X2\x^T = —(1 + RX2 + P-x\)x2 + 
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Figure 7: The red arcs are the trajectory of the system (|5.6()|1 with initial value (0,1.2) and the blue arcs 
are the orbit with initial value (0,1.5). Through simulation, we observe that the trajectories approach 
to the periodic solution of (|5.60p as time increases. 


where R = 0.9 and T = {x|xi = 0, X 2 < 0}. It is easy to see that system ^5.69\) is of the form ^5.57\) and 
^{xi^X 2 ) = Xi = 0. The system has a periodic solution 


T^(i) = (1 -h (1 + ^) cos(t), -(1 + fi) sin(i)), 


(5.70) 


where i G M for p G (—2,0]. 

The generating system of i5.69\) has the following form 


(5.71) 


Xi = X2, 

x '2 = —0.0001[x2 + {xi — 1)^ — 1]X2 — Xi + 1, 

Aa;2|xGr = -(1 + Rx2)x2, 

and admits the periodic solution = (l+cos(t), — sin(t)). By means of the equality (V$(a:*), f{x*)) = 

((1,0), (0,1)) = 0 with X* = (0,0) € 9r, it is easy to say that x* is a grazing point o/'I'o(i). 

Let us start with the linearization of system (I5.71|) around the periodic solution 'I'o(0- Consider a 
near solution y[t) = y(t, 0, y* + Ay), where Ay = (Aj/i, Ay 2 ), to the periodic solution 4'o(i). Assume that 
y{t) satisfies condition (iVl), and it meets the surface of discontinuity T at the moment t = and at the 
point y = y{£^, 0, y* + Ay). Considering the formula (13.101) for the transversal point y = {yi, 2 / 2 ), the first 


component 


dT-jy) 

dy^ 


can be evaluated as 


drjy) 

dyi 


= — — . From the last equality, the singularity is seen at the 


y2 


grazing point. By taking into account (13.171) with (15.711) and 


drjy) 

dyf ’ 


we obtain that 
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(5.72) 


dWM 

dy° 


Ry2 -1 

-0.0001i?(yi + (yi - 1)2 - 1) - 2i?(0.0001(y2 + - 1)2 - 1)) 


Similarly, taking into account the formula (I3.22|l . one can evaluate that 
(13.2311 imply 


dT{y) 

dy°2 


0. This and formula 


dW^{y) 

dy°2 


0 

-2i?y2 


Joining (15.721) and (15.7311 . the matrix Wiy{y) can he obtained as 


(5.73) 


Wiy{y) 


Ry2 - 1 0 

-0.0001i?(yi + (yi - l)^ - 1) - 2i?(0.0001(yi + (yi - 1)^ - 1)) -2Ry2 


(5.74) 


The last expression implies continuity of the partial derivatives near the grazing point. This validates 
eondition {A2). 

Then, evaluating the matrix in (15.741) at y = y* = (0,0), it is easy to obtain 


W.yiy*) 


-1 0 
0.0003i? 0 


(5.75) 


and 


O 


2 , 




O.OOOSi? 0 


if (N^) is valid, 

if (.^2) is valid. 


To verify condition (A3), let us specify the region 


H = {(,yi,y 2 )\y 2 < - {yi - l)^, 0 < yi < 1}. 


(5.76) 


For the grazing solutionwe have that6i.\.i—6i = 27r. Consider a near solution y{t) = (yi(t),y 2 {t)) = 
y{t,0,y) to 'l'(i). To satisfy the condition (Nl), take y = {yi,y 2 ) S H. The orbit of y(t) is below the 
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grazing orbit. Fix points y = {yi,y 2 ) G H and ^ = ('i/'i, ^” 2 ) of the orbits y{t) and 'l'o(i), respectively such 
that 0 < t/i = '01 < 1 and ip 2 < 0. Since of the equation y[ = 1/2, the speed of yi{t) at (?/i,t/2) is larger 
than the speed o/4'i(t) at (0i,02)- Consequently, one can find that T{y) < j < 27r for y £ H. Thus, the 
condition (^3) is valid and Lemma \3.1\ verifies the condition (C). 

It is easy to demonstrate that the condition {AT) is valid such that near solutions to the grazing 
one are either continuous or discontinuous. That is, they don’t intersect the line of discontinuity T or 
intersect it permanently near to the grazing point and by means of the formula (15.761) . the linearization 
system for (15.711) around the grazing cycle 'I'o(0 consists of the following two subsystems 


u'l = U2, 

u '2 = —0.0001 sin(2t)Mi + 0.0002 sin^(t)M2, 


(5.77) 


and 


Ui = U2, 

U 2 = —0.0001 sin(2t)ui + 0.0002 sin^(t)u 2 , 

-1 0 
0.00037? 0 


^^|27r2 — 


U. 


(5.78) 


The system (15.771) + (15.781) is {2 tt, 1) periodic. The Floquet multipliers of system (15.771) + (15.781) are 
= 1, = 0.939, = 1, p^^ = 0.912. Thus, condition (^15) is validated. Moreover, the conditions 

(Cl) — (C7) and (Al), {A2) can be verified utilizing similar way presented in Examvle \S. 1[ Consequently, 
Theorem \4.1\ authenticates that the grazing periodic solution (cycle), To{t) of the system (15.711) is orbitally 
stable. The simulation results demonstrating the orbital stability ofT^ft) are depicted in Figure\^ 



Figure 8: The grazing cycle of system (|5.71l) is in red. The blue arcs are the trajectory of the system 
with initial point (0.5,1.2) and the green continuous orbit is with initial value (0.1,0). They demonstrate 
stability of the grazing solution. 


Next, we will investigate two sorts of periodic solutions of system Ii5.69\) with a period AT near to 2 tt. 
The first one is continuous and the second admits discontinuities once on a period. For those solutions, 
corresponding linearization systems around the grazing cycle To(t) are (I5.77|) and (15.781) . respectively. 
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Let us start with the continuous periodic solutions of i5.69\) . For continuous periodic solution, we will 
consider the linearization system (15.771) . 

To apply Theorem \5.1{ denote 'I'o(O) = (Ci,0). That is, consider (2 = 0. Then, applying the above 
discussion, obtain that the Poincare condition admits the form of the following equations, 


-xi=0, 


(5.79) 


^2(^, Cl, M) = X2{.3^, Cl, M) = 0. 

Because solutions of the system (15.711) have continuous derivatives with respect to the time, phase 
variables and parameters, we can calculate the following determinant 


9^i(u;,Ci°,0) Ci°,0) 

d3T dx\ 

9 ^ 2 ( 0 ;, Ci°,o) Ci°,o) ■ 

dx\ 

First, we need the monodromy matrix of the system (15.771) . It is 


(5.80) 


0.939 -0.0001407 

-0.0003165 1 


(5.81) 


It is easy to see that first column of the determinant i5.80\} is computed by utilizing h5.71\) and the 
second column is evaluated by means of the first column of the matrix i5.81\) . From this discussion, one 
can obtain that the determinant (15.801) is equal to 


0 -0.061 
1 -0.0003165 


0.061 7 ^ 0 . 


(5.82) 


Thus, in the light of Theorem \5.1[ we can conclude that for sufficiently small |/i| there exists a unique 
periodic solution of the system 


x'l = X2, 

x'2 = —0.0001[x2 + (xi — 1 )^ — (1 + h)'^]x2 — Xi + 1 . 


(5.83) 


It is exactly the cycle h5. 701) with a period S' = 2?!. If p, < 0, the solution is separated from the set 
r. Consequently, it is a periodic continuous solution of the equation i5.69\) . It is orbitally stable by the 
theorem for continuous dynamics since of the continuous dependence of multipliers on the parameter. 
The function "i!^(t),fi > 0, intersects T and can not be a solution of equation Thus, the system 

does not admit a continuous periodic solution near to 'I'o( 0 , */ parameter is positive. 

Considering those solutions which have one moment of discontinuity in a period, one can find that 


38 































the corresponding linearization of'^o{t) is the system (15.781) . 
The monodromy matrix of (15.781) can be evaluated as 


0.939 -0.00052 

-0.000427 1 


(5.84) 


It can be easily observed that the discontinuous solution intersects the line xi = 0. For this reason, one 
can specify the first coordinate of the initial value as fi = (^i = 0. In the light of these discussions and 
the formula (I5.58|) . the following equations are obtained: 


^i(^,0,C2,/i)=a;i(=(^,0,C2,M)=0, 

^2(^, 0, C2, /i) = 0, C2, M) - C2 = 0. 


(5.85) 


Then, taking the derivative of the system (15.851) with respect to and (^ 2 , o-nd calculating it at 7 = ui, 
(2 = C 2 — 0, and for = 0, the following determinant is obtained 


a^i(c^,o,C2°,o) 

5^i(w,0,C2°,0) 



d7 

3(2 


O -0.0006 

a^2(c^,0,C2°,0) 

5^2(w,0,C2°,0) 


1 0.0009 

d7 

d(2 




(5.86) 


Thus, condition (A6) holds. Then, utilizing Theorem \5.1[ it is easy to conclude that for sufficiently small 
p, there exists a unique periodic solution of the system S5.69\) with a period « 2tt. It is true that for 
positive as well as negative p. Moreover, these solutions are orbitally asymptotically stable because of the 
continuous dependence of solutions on parameter and initial values and they meet the discontinuity line 
transversally. 

For each fixed p 0, solutions near to the periodic ones intersect the line of discontinuity F transver¬ 
sally once during the time approximately equal to the period. That is, the smoothness which is requested 
for the application of the Poincare condition is valid, since the smoothness for the grazing point has 
already been verified. It is clear that there can not be another solutions with period close to 2 tt. Thus, 
one can make the following conclusion. The original system i5.69\} admits two orbitally stable periodic 
solutions, continuous and discontinuous, if p < 0. There is a single orbitally stable continuous solution 
(grazing) if p = 0. Additionally, there is a unique discontinuous orbitally stable periodic solution for 
positive values of the parameter. Consequently, grazing bifurcation of cycles appears for the system with 
small parameter. 

We have obtained regular behavior in dynamics near grazing orbits by the Poincare small parameter 
analysis. Nevertheless, outside the attractors irregular phenomena may be observed. 

In Figure\^ the solutions of the system h5.69\) with parameter p = —0.2 are depicted through simula¬ 
tions. The red arcs are the trajectory of the system i5.69\) with initial value (0.7,0.05) and the blue arcs 
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are the trajectory of the system iS.Offl) with initial value (0.4,0.05). It is seen that both red and blue tra¬ 
jectories approach the discontinuous periodic solution of i5.69\} . as time increases. So, the discontinuous 
cycle is orbitally stable trajectory. Moreover, the green one is a continuous periodic trajectory of Ii5.69\) 
with initial value (0, 0.05) and it is orbitally asymptotically stable. To sum up, there exists two periodic 
solutions of \5.69\} for the parameter p, = —0.2, one is continuous, the other one is discontinuous and 
both solutions are orbitally asymptotically stable. 



Figure 9: The blue, red and green arcs constitute the trajectories of system (|5.69H with p = —0.2. The 
first two approach as time increases to the discontinuous limit cycle and the third one is the continuous 
limit cycle itself. 

In Fig. \1(A the red arcs are the orbit of the system with initial value (0, 0.1) and the blue arcs are the 
trajectory of it with initial value (0,0.4). Both trajectories approach to the discontinuous cycle of system 
h5.69\) . as time increases. Thus, Fig. \1(A illustrates the existence of the orbitally stable discontinuous 
periodic solution if p = 0.2. 



Figure 10: The red and blue arcs constitute the trajectories of the system (15.6911 with p = 0.2. Both 
orbits approach to the discontinuous limit cycle, as time increases. 


6 Conclusion 

In literature, the dynamics in the neighborhood of the grazing points d-nD], [13], [ig, eg, ESj-isg 
is generally analyzed through maps of the Poincare type. The main analysis is conducted on complex 
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dynamics behavior such as chaos and bifurcation Hen], na, im, na, issi-Eni. However, there is still 
no sufficient conditions for the discontinuous motion to admit main features of dynamical systems : the 
group property, continuous and differentiable dependence on initial data and continuation of motions, 
which are useful for both local and global analysis. Variational systems for grazing solutions have not 
been considered in general as well as orbital stability theorem and regular perturbation theory around 
cycles, despite, particular cases can be found in specialized papers. See, for example, m- To investigate 
these problems in the present paper, we have applied the method of B— equivalence and results on 
discontinuous dynamics developed and summarized in [1]. In our analysis the grazing singularity is 
observed through the gradient of the time function t(x), since some of its coordinates are infinite. We 
have found the components of the discontinuous dynamical system that is the vector field, surfaces of 
discontinuity and the equations of jump such that interacting they neutralize the effect of singularity. 
Then, we linearize the system at the grazing moments and this brings the dynamics to regular analysis and 
make suitable for the application. By means of the linearization, the theory can be understood as a part 
of the general theory of discontinuous dynamical system. Thus, we have considered grazing phenomena 
as a subject of the general theory of discontinuous dynamical systems [T], discovered a partition of set 
of solutions near grazing solution such that we determine linearization around a grazing solution is a 
collection of several linear impulsive systems with fixed moments of impulses. This constitutes the main 
novelty of the present paper. To linearize a solution around the grazing one, a system from the collection 
is to be utilized. This result has been applied to prove the orbital stability theorem. The way of analysis 
in m-m continues in the present paper and it admits all attributes which are proper for continuous 
dynamics |23| . That is why, we believe that the method can be extended for introduction and research 
of graziness in other types of dynamics such as partial and functional differential equations and others. 
Next, we plan to apply the present results and the method of investigation for problems initiated in 

H-nni, m-m, m- 
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